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^ , Abstract 

fN| ' Primordial nucleosynthesis as well as anisotropies in the cosmic mi- 

crowave background radiation imply that the total amount of baryons 
in the Universe largely exceeds the visible contribution, thereby mak- 
ing a strong case for baryonic dark matter. Moreover, certain recent 
developments lead to a consistent picture of the dark baryon budget 
in the present-day Universe. Accordingly, dark baryons are mostly 
locked up in galactic halos - which are anyway dominated by non- 
baryonic dark matter - and a sizable fraction of them consists of gas 
clouds. While a priori various forms of baryonic dark matter in galax- 
ies can be conceived, observational constraints rule out most of the 
possibilities, leaving brown dwarfs and cold gas clouds mostly made of 
2 1 H2 as the only viable candidates (besides supermassive black holes). 

c/2 , So, it looks natural to suppose that baryonic dark matter in galaxies 

is accounted for by dark clusters made of brown dwarfs and cold H2 
clouds. A few years ago, it was shown that indeed these dark clusters 
are predicted to populate the outer halos of normal spiral galaxies 
by the Fall-Rees theory for the formation of globular clusters, which 
was based on the standard cold dark matter paradigm described in 
Blumenthal et al. 1984 Nature 311, 517. We review the dark clus- 
ter formation mechanism, and argue that its qualitative features are 
expected to remain true even in the contemporary picture of galaxy 
formation. We also discuss various ramifications of the dark cluster 
scenario in question, paying particular attention to its observational 
implications. One of them - the diffuse gamma-ray emission from the 
Milky Way halo - appears to have been confirmed by the discovery of 
Dixon et al. 1998 New Astronomy 3, 539. Whether this is actually 
fact or fiction only the future satellite missions AGILE and GLAST 
will tell. 
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1 Introduction and motivation 



All available observational evidence - including flat rotation curves of spi- 
ral galaxies, X-ray emission from elliptical galaxies and clusters of galaxies, 
gravitational lensing by galaxy clusters, velocity fields in deep galaxy sur- 
veys, light curves of high-redshift Type-la supernovae and cosmic microwave 
background radiation (CMBR) anisotropies - invariably lead to the conclu- 
sion that the Universe is dominated by dark matter. Denoting (as usual) 
by Q the actual density of the present-day Universe in units of the critical 
density, the situation can be summarized as follows. 

Measurements of anisotropies in the CMBR yield both the total value 
f2 ~ 1 from the location of the first acoustic peak 1 (De Bernardis et al. 
2002; Balbi et al. 2000; Pryke et al. 2002) and the baryonic contribution 
from the height of the second peak 2 Qj, ~ 0.04 (Sievers et al. 2002; Stompor 
et al. 2001). Quite remarkably, the latter value is in good agreement with 
two independent estimates. One arises by combining measurements of the 
primeval abundance of light elements (especially deuterium) with theoretical 
predictions from primordial nucleosynthesis (Buries et al. 2001; O'Meara 
et al. 2001). The other follows from the observed features of high-redshift 
Lyman-a forest absorption lines of neutral hydrogen (Rauch et al. 1998) 
observed in the spectra of background quasars 3 . However, this is not the 
end of the story, since estimates of the total matter distribution - based on 
the spectral features of the CMBR, on deep galaxy surveys and on various 
studies of clusters of galaxies - yield Vl m ~ 0.3 (Turner 2002). Evidently, 
matter alone - as clustered on various cosmic scales - fails to account for the 
total energy budget of the Universe, leaving a gap of Q — Q m ~ 0.7. Actu- 
ally, the situation is even more puzzling, since observations of high-redshift 
Type-la supernovae entail that the present cosmic expansion is accelerated. 
Nevertheless, this behaviour is explained by a nontrivial vacuum (unclustered 
energy) with negative pressure and — 0.7 4 (Riess et al. 1998; Perlmutter 

1 Observe that this result agrees with the natural implication of cosmic inflationary 
scenarios (Liddle & Lyth 2000). 

2 Behind this second conclusion there is the assumption that primordial fluctuations are 
adiabatic, which is just what cosmic inflation entails. The value Hq ~ 70 km s _1 Mpc^ 1 
for the Hubble constant is adopted throughout. 

3 Lyman-a forest lines will be addressed in more detail in Sect. 3. 
4 This vacuum is described by a cosmological constant A. 
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et al. 1999). Notice that such a vacuum provides just the missing contri- 
bution Q — Q m ~ 0.7 mentioned above, and so the overall picture is indeed 
remarkably consistent. 

The foregoing discussion implies Vt b /Vt m ~ 0.13, thereby showing that 
most of the dark matter is nonbaryonic. This is in agreement with contem- 
porary theories of galaxy formation, which demand that the observed dark 
galactic halos should be dominated by nonbaryonic matter. 

Still - as we will see below - half of the total amount of baryons in the 
present-day Universe remains observationally unaccounted. Hence, there is 
also a baryonic dark matter problem, which is the main topic addressed in 
this paper. 

An inventory of the baryonic content of the various cosmic structures has 
been attempted by Persic and Salucci (1992), Gnedin and Ostriker (1992), 
Bristow and Phillipps (1994), and Fukugita, Hogan and Peebles (1998) by 
combining available observational data with the whole body of theoretical 
knowledge. According to the latter authors, the baryon budget in the present- 
day Universe yields Q' b ~ 0.02 and is dominated by hot gas in groups and 
clusters of galaxies, while the equivalent cosmological contribution from lu- 
minous baryons 5 in galaxies is as small as Q! bg ~ 0.004. 

Yet, we know that the observed CMBR anisotropies, primordial nucle- 
osynthesis and high-redshift Lyman-a forest absorption require Qb — 0.04. 
Where are the dark baryons"! 

An answer to this question is provided by hydrodynamic simulations of 
currently favoured structure formation models (Cen & Ostriker 1999; Dave 
et al. 1999). They indicate that - in the present-day Universe - baryons 
are divided almost equally among three phases: (i) warm gas (T = 10 5 — 
10 7 K), forming transient filamentary structures in intergalactic space; (ii) 
cool gas (T < 10 5 i^), smoothly diffuse in intergalactic space; and (iii) cold 
gas and stars associated with galaxies (Dave et al. 2001). Although the 
uncertainties in the simulations prevent a precise estimate of the fractional 
abundance of baryons in the three phases, the equivalent cosmological density 
in each one is roughly 0.013. This means that the equivalent cosmological 
contribution from galactic baryons should be Q% ~ 0.013. Because we have 
seen that luminous baryons in galaxies only give Q' b ~ 0.004, we are led to 

5 Luminous baryons include stars, atomic and molecular hydrogen whose existence is 
inferred from detection of the emitted electromagnetic radiation at any wavelenght. 
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the conclusion that in the low-redshift Universe most of the galactic baryons 
are dark. 

Actually, a strong case for a somewhat larger amount of baryonic dark 
matter in galaxies has recently been made by Kochanek (2001). He has 
investigated the mass function of the dark matter halos - as predicted by 
the current cold dark matter paradigm - by using certain dynamical probes 
(Kochanek & White 2001). More specifically, Kochanek has shown that 
the distribution of gravitational lens separations and the circular velocity 
function predicted by semi-analytic models of hierarchical clustering are both 
in agreement with observations only provided that Qbg — 0.02 6 . Basically 
the same conclusion has been reached by Klypin, Zhao and Somerville (2001) 
via N-body simulation modelling of the Milky Way and Andromeda galaxy. 
These authors found that the observed properties are reproduced only if 
50 % — 75 % of the galactic baryons are nonluminous. Consequently, galactic 
halos are the main repository of dark baryons at low redshift 7 8 . 

As a matter of fact, this circumstance is in remarkable agreement with 
the claim of Lanzetta and collaborators (Chen et al. 1998, 2001) that low- 
redshift Lyman-a forest absorption lines of neutral hydrogen are produced by 
baryonic gas clouds located inside the dark galactic halos. Notice that these 
gaseous baryons are missing in the above-mentioned inventories. 

Apparently, a consistent picture of the dark baryon budget in the low- 
redshift Universe emerges from the previous considerations. Dark baryons 
are mostly locked up in galactic halos, and a sizable fraction of them consists 

6 We remark that no serious conflict exists between this conclusion and the outcome of 
hydrodynamic simulations, owing to their intrinsic uncertainties. 

7 Of course, this does not undermine the fact that galactic halos are dominated by 
nonbaryonic dark matter. 

8 Strictly speaking, there is a logical gap here, which should be discussed. Contrary 
to hydrodynamical simulations, N-body simulations and semi-analytical models of galaxy 
formation neglect feedback processes, which give rise to galactic outflows. Therefore, the 
above results of Kochanek and Klypin, Zhao and Somerville do not prevent ejection of a 
large amount of baryons during the galaxy evolution. Nevertheless, there are good reasons 
to believe that most of the baryons should still lie inside galactic halos today. For, galac- 
tic winds powered by supernova explosions appear unable to drive a huge amount of gas 
out of galactic halos (Mac Low & Ferrara 1999). In addition, a large amount of ejected 
baryons from spiral galaxies would produce an equivalent cosmological density of the in- 
tergalactic medium larger than the one predicted by the above-mentioned hydrodynamical 
simulations. 
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of gas clouds. 

Several kinds of astronomical objects - like brown dwarfs 9 , red dwarfs, 
white dwarfs, neutron stars, black holes and gas clouds - may account for 
baryonic dark matter in galaxies (Carr 1994), but observational constrains 
rule out most of them. Before a star becomes a white dwarf a neutron star or 
a black hole, it goes through phases in which heavy elements are produced and 
ejected into interstellar space and its luminosity is greatly enhanced. Existing 
upper limits on the metallicity of the environment and on the background 
light then rule out these objects as candidates for baryonic dark matter 10 
(Carr 2000; Madau et al. 2000). A similar conclusion holds for red dwarfs 
as well, because they would give rise to infrared fluxes stronger than are 
detected (more about this, in Sect. 3). As far as gas clouds are concerned, 
only cold gas mostly made of H 2 is a viable possibility, for otherwise an 
unseen excess of photons at some wavelenght (21 cm for neutral hydrogen) 
would be produced 11 . So - with the exception of supermassive black holes 
- only brown dwarfs and cold H 2 clouds turn out to be realistic forms of 

9 We recall that brown dwarfs are stellar-like objects in which ordinary nuclear reactions 
do not occur in the core, because their central temperature never exceeds the hydrogen- 
burning threshold (for a review, see Kulkarni 1997; Basri 2000; Chabrier & Baraffc 2000). 
In the case of solar metallicity, the brown dwarf mass range is 0.01 — 0.08 M , but for 
primordial metallicity the upper limit becomes M < 0.11 M Q (D'Antona 1987). For a 
long time, brown dwarfs were regarded merely as a theoretical possibility, until they were 
first discovered in 1995 (Rebolo et al. 1995; Nakajima et al. 1995). Further observations 
have shown that brown dwarfs form like ordinary stars. Similarly to what happens in 
very low-mass stars (Drake et al. 1996), it has been argued that also in brown dwarfs a 
turbolent dynamo process should heat the corona, which consequently emits soft X-rays 
(Kashyap et al. 1994; Neuhauser et al. 1999). This point will be important in some later 
considerations. 

10 The metallicity constraint is evaded by black holes with mass > 200 M , because they 
would ingurgitate their whole parent stars. However, only supermassive black holes - with 
mass > 10 5 M Q - turn out to be viable dark matter candidates, for otherwise the parent 
stars would violate background light constraints. 

11 As H2 does not possess a permanent electric dipole moment, the first transition line to 
be detected in emission occurs by quadrupolc radiation at 512 K above the ground state. 
Thus, cold H 2 clouds do not emit any radiation (Combes & Pfenniger 1997). Notice that 
astronomers tend to infer the existence of H2 by using CO as a tracer. Needless to say, 
this practice makes sense only when the H2/CO conversion ratio is under control (like in 
the Milky Way disk), but this is not the case in galactic halos. The possibility that a large 
amount of H 2 may be present in the Universe was first suggested by Zwicky (1959). 
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baryonic dark matter in galaxies 12 13 14 . 

Our aim is to discuss the possibility that the dark galactic baryons indeed 
consist of dark clusters - made of brown dwarfs and cold H 2 clouds - lurking 
in the halos of normal spiral galaxies. 

More specifically, we will focus our attention on a model of this kind, 
which was developed from first principles by De Paolis, Ingrosso, Jetzer and 
Roncadelli (1995a,b, 1996, 1998a) within the standard cold dark matter pic- 
ture of galaxy formation. We stress that an almost identical scenario has 
been independently proposed and investigated by Gerhard and Silk (1996). 
Somewhat similar ideas have been put forward by Ashman and Carr (1988), 
Ashman (1990), and Fabian and Nulsen (1994, 1997). Slightly different bary- 
onic pictures have been suggested by Pfenniger, Combes and Martinet (1994), 
Wassermann and Salpeter (1994), Gibson and Schild (1999a,b), Sciama (2000), 
and Lawrence (2001). 

Here, we offer an updated review of the above-mentioned dark matter 
model. 

The paper is organized as follows. In Sect. 2, we discuss the formation 
mechanism (Subsect. 2.1) and the main properties (Subsect. 2.2) of the dark 
clusters. Next, we address some of their observational implications in Sect. 
3. Finally, we summarize our conclusions in Sect. 4. Suggestions for further 
work will be scattered throughout the paper. 

12 Of course, a small fraction of red dwarfs, white dwarfs, neutron stars and black holes 
can be present in the galactic halos, but they cannot account for the bulk of baryonic dark 
matter. 

13 A strong constraint on the brown dwarf fraction in galactic halos - which would prevent 
them from playing any role for dark matter - was derived by Graff and Freese (1996a) 
under the assumption that the initial mass function is spatially homogeneous. However, 
within the dark matter model to be discussed below the initial mass function depends on 
the galactocentric distance, and brown dwarfs can copiously form only at galactocentric 
distances larger than 10 — 20 kpc. Hence, the conclusion of Graff and Freese (1996a) is not 
valid in the present context. 

14 We anticipate (from Sect. 3) that old brown dwarfs are not ruled out by current 
infrared searches (especially if they are clustered). 
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2 Dark clusters in galactic halos 



In the past, brown dwarfs have been repeatedly recognized as an important 
contributor to baryonic dark matter in galaxies (Carr 1994). Much like or- 
dinary stars, brown dwarfs are expected to form in clusters, thereby giving 
rise to dark clusters (Carr & Lacey 1987; Silk 1991; Rix & Lake 1993; Kerins 
& Carr 1994; Moore & Silk 1995). However, insufficient attention was paid 
to the fact that - even in this case - the star-formation mechanism is highly 
inefficient, in that most of the original gas fails to get transformed into stars, 
thereby remaining trapped inside the dark clusters. As we shall see, this gas 
plays a crucial role, since it can both change the properties of the brown 
dwarfs and make the dark clusters directly visible in the gamma-ray band 
and indirectly observable along the line of sight to a distant quasar. 

2.1 Formation of the dark clusters 

Starting point of the considered dark matter model was the Fall-Rees (1985) 
theory for the formation of globular clusters in galactic halos. Indeed, it 
was shown (De Paolis et al. 1995a,b, 1996, 1998a) that the Fall-Rees theory 
automatically predicts - without any further physical assumption - that dark 
clusters made of brown dwarfs and cold H 2 clouds should form in the outer 
halos of normal spiral galaxies, namely at galactocentric distances larger than 
10 - 20 kpc. 

The Fall-Rees theory was formutated within the standard cold dark mat- 
ter picture of galaxy formation described in Blumenthal et al. (1984). How- 
ever, in the last few years numerical simulations of structure formation have 
led to a considerable improvement in the understanding of hierarchical clus- 
tering, with mergers playing an increasing role. As a result, some features 
of the above-mentioned cold dark matter paradigm turned out to be incor- 
rect and had to be revised. Therefore, it is not clear whether the present 
dark matter model remains true. Although this is at the moment an open 
question, we will argue - on the basis of the following discussion - that the 
qualitative features of the dark matter scenario under consideration are likely 
to remain unchanged even in contemporary models of galaxy formation. 

Consider a protogalactic cloud with mass M ~ 1O 12 M , made of both 
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nonbaryonic and baryonic matter 15 . We imagine that its inner part eventu- 
ally produces - upon gravitational collapse - the luminous component of a 
normal spiral galaxy whereas its outer region gives rise to the corresponding 
dark halo. 

The Fall-Rees theory rests on the assumption that the outer part of the 
protogalactic cloud is modelled by a singular isothermal sphere (SIS), so that 
the observed flat rotation curve is recovered. Accordingly, the overall density 
profile is p (r) = V max / 'AnGr 2 , where V max is the value at which the rotation 
velocity becomes constant at a large galactocentric distance r. Typically one 
finds Vmax ~ 200 km s _1 , and we will choose V max = 220 km s _1 , which is the 
value appropriate to the Milky Way. Hence, we get 

Po (r)^ 0.9 (^] M pc- 3 . (1) 



Correspondingly, the free-fall time tff ~ 0.5 a 1 l 2 (Gp ) x l 2 becomes 

%^7.4.10 6 a-V2 (j-^j yr , (2) 

where a is a constant of order unity which reflects the actual mass distribu- 
tion (a pure SIS yields a — 1). Moreover, the mean square velocity of any 
baryon in the considered region is v 2 = 3V max /2 ~ 7.3 • 10 4 km 2 s~ 2 . There- 
fore the temperature of the diffuse baryonic gas 16 is T d ~ 1.7 • 10 6 K. One 
can easily see that this temperature is quite close to the virial temperature, 
and so the diffuse gas is initially in virial equilibrium. 

Since the diffuse gas tends to cool, in order to understand its further 
behaviour one has to compare the cooling time t coo i with the free-fall time 
tff as given by eq. (2) 17 . One can show that at the above temperature both 
Bremsstrahlung and ion-recombination processes contribute to the cooling 
rate of the diffuse gas, and the resulting cooling time is 

t cool ~ 5.5 • 10- 18 (^jf] yr , (3) 



15 Wc suppose that the baryonic matter has primordial composition. 

16 As we will be explicitly concerned with the baryonic gas component only, the attribute 
baryonic will be implicitly understood. Also, the meaning of the attribute diffuse will 
become clear later on. 

17 We recall that a cloud collapses in free-fall and fragments when t coo i < tff, whereas 
the collapse proceeds quasi-statically for t coo i > tff. 
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where pd denotes the density of the diffuse gas. The comparison between eqs. 
(2) and (3) entails 

tcool P* 



tff Pd 



(4) 



where we have set 

p* 



1.1 - lO- 2 ^ 1 / 2 (^] M Q pc' 3 . (5) 



On account of eq. (4), we would conclude that the diffuse gas undergoes 
quasi-static collapse as long as its density is sufficiently small - namely for 
Pd < P* - since then t coo i > tff. Still, as soon as pd > p* the collapse enters 
the free-fall regime. Actually, in the latter case we would expect a monotonic 
increase of pd, and - thanks to eq. (4) - a steady decrease of the t coo i/tff 
ratio, which would signal efficient fragmentation and star formation within 
the whole protogalactic cloud. 

But in reality the situation is more complex, since the above picture 
would be correct only within a perfectly homogeneous medium. In fact, Fall 
and Rees have shown that the unavoidable presence of density fluctuations 
during the collapse can totally upset such an expectation. Basically, this is 
due to the fact that density fluctuations produce a thermal instability within 
the diffuse gas. To see this, consider eq. (3): an overdense region cools more 
rapidly than average, but then compression from the surrounding hotter gas 
leads to a further increase of the density. Thus, density fluctuations give 
rise to a substructure within the protogalaxy, made of denser and cooler gas 
bubbles embedded in the diffuse gas. Besides, Fall and Rees have shown that 
the bubble condensation implies the condition 

h^L ~ i ( 6 ) 

to hold true for the diffuse gas. This fact is extremely important, for three 
different reasons. First, the diffuse gas actually stays in virial equilibrium at 
T d ~ 1.7 • 10 6 K. Second, there is no star formation within the diffuse gas. 
Third, eqs. (4), (5) and (6) imply that the density profile of the diffuse gas 
is 

AfcCrO-l.l.lO-W 2 (^) M pc" 3 . (7) 
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The bubbles - which are pressure- confined by the diffuse gas - are origi- 
nally made of ionized hydrogen and helium (plus their respective electrons). 
Assuming that this plasma is in thermal equilibrium, the cooling brought 
about Bremsstrahlung and ion-recombination processes is operative only at 
temperatures larger than ~ 10 4 K, since at lower temperatures the gas be- 
comes neutral and the corresponding cooling rate suddenly drops to zero. 
Hence, we expect that inside the bubbles hydrostatic equilibrium sets in when 
the temperature reaches ~ 10 4 K. In such a situation, the baryonic compo- 
nent of the protogalaxy is made of bubbles with T b ~ 10 4 K embedded in the 
diffuse gas at Td ~ 10 6 K. 

Let us denote by p&(r) the (constant) density inside a bubble located at 
galactocentric distance r. At hydrostatic equilibrium, the pressure p b just 
inside the considered bubble coincides with the pressure pd just outside it 
(within the diffuse gas). This condition entails 




where ~ 10~ 24 g and mb ~ 2 • 10~ 24 g are the mean particle mass of a ion- 
ized and neutral gas having primordial composition, respectively. Combining 
eqs. (7) and (8) together, we find 

p b (r) ~ 3.6 a 1 / 2 M pc" 3 . (9) 

Therefore the Jeans mass and the Jeans radius for the bubble in question are 

and 

respectively, where k B denotes the Boltzmann constant. 

However, the above assumption - namely thermal equilibrium of the 
plasma inside the bubbles - turns out to be violated, since cooling occurs 



10 



very rapidly (Kang et al. 1990). The resulting out- of- equilibrium ion re- 
combination entails a substantial ionization fraction at Tj, < 10 A K (which 
would otherwise vanish at thermal equilibrium). This circumstance is in 
fact irrelevant for T& > 10 4 K, but alters dramatically the above conclusions 
for T5 < 10 4 K. Indeed, the presence of a nonnegligible amount of ions at 
T5 < 10 A K gives rise to the formation of molecular hydrogen via the reac- 
tions 18 H + p->H£ + <y,H + e^H- + ~f and H 2 + + H -> H 2 + p , 
H + H" — > H 2 + e (notice that ions only act as catalysts). 

The presence of H 2 in the bubbles brings about a further cooling, owing 
to photon emission in roto-vibrational molecular transitions. We stress that 
this process is very efficient and can cause a temperature drop down to 
T b ~ 10 K. Evidently in such a situation t coo i <C £//, and so the bubbles 
collapse in free-fall and fragment. When the number density in the bubbles 
exceeds ~ 10 8 cm~ 3 - which corresponds to pi> ~ 2.8 • 10 6 M & pc~ 3 - a further 
H 2 production takes place via the three-body reactions H + H + H — > H 2 + H 
and H + H + H 2 — > H 2 + H 2 , as realized by Palla, Salpeter and Stahler (1983). 
At variance with the former reactions, the latter ones do not involve ions as 
catalysts. As a consequence, the bubbles get transformed almost entirely into 
H 2 . Besides modifying their internal composition, this circumstance makes 
cooling even more efficient. 

Still, a priori nothing ensures that H 2 - once produced - will survive: 
because of its fragility, its actual existence crucially depends on the environ- 
mental conditions. 

In fact, in the central region of the protogalaxy an AGN (Active Galactic 
Nucleus) along with a first population of massive stars (population III) 
are expected to form. They act as strong sources of ultraviolet radiation, 
which dissociates the H 2 molecules up to a critical galactocentric distance r* 
(because the radiative flux decreases as the galactocentric distance increases). 
It is not difficult to estimates that the H 2 destruction should occur for r < r*, 
with = 10 — 20 kpc. Consequently, the further evolution of the inner 
region of the protogalaxy (r < 10 — 20 kpc) will be totally different from that 
of the outer part (r > 10 — 20 kpc). Actually, since the spheroid of radius 

18 As already pointed out, H2 has no permanent electric dipole moment, and so the 
standard reaction H+H — > H 2 +j can only occur on dust grains. But since the composition 
of the protogalaxy is supposed to be primordial, no grains exist. Hence, the considered 
reaction is not operative. We stress that this fact implies that CO is not a tracer of H 2 in 
the present context. 
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10 — 20 kpc gives rise to the inner halo of a normal spiral galaxy, we will 
distinguish between inner and outer halo. 

Owing to the above considerations, we are now in position to summarize 
the dynamics of the halo as follows. 

• Inner halo (Fall & Rees 1985) - The lack of H 2 prevents the bub- 
bles from cooling down to temperatures smaller than ~ 10 4 K, and so 
they stay for a long time in hydrostatic equilibrium, which coincides 
with virial equilibrium (since t coo i tff for T b < 10 4 K). As a result, 
their mass and radius are given (respectively) by the corresponding 
Jeans values (10) and (11). Indeed, the main result of the Fall-Rees 
theory is that these values are in agreement with those observed for 
globular clusters provided that mass loss during formation and evolu- 
tion is taken into account. Subsequently the ultraviolet flux decreases, 
thereby allowing for the formation and survival of H2. Accordingly, the 
bubbles can further cool, collapse in free-fall and fragment, ultimately 
producing ordinary stars. Thus, the Fall-Rees theory provides a nat- 
ural explanation for the formation of globular clusters in the halos of 
normal spiral galaxies. 

• Outer halo (De Paolis et al. 1995a,b) - The presence of H 2 allows the 
bubbles to cool down to temperatures much lower than 10 4 i\", thereby 
collapsing in free-fall and fragmenting. As is well known, the collapse 
becomes adiabatic - and eventually stops - when a fragment becomes 
optically thick. Palla, Salpeter and Stahler (1983) have shown that this 
occurs when the fragment Jeans mass is as low as 10~ 2 — 10 _1 M e . So, 
clusters of brown dwarfs in the mass range 10~ 2 — 10 _1 M should form 
in the outer halos of normal spiral galaxies. 

Notice that the present dark matter model gives rise to an initial mass func- 
tion in the halo that depends on the galactocentric distance 19 20 . 

19 Evidence for a spatially varying initial mass function in the Milky Way disk has been 
reported by Taylor (1998). 

20 Because of this fact, no extrapolation from the initial mass function in the disk can 
yield informations about the initial mass function in the halo. Consequently, the bound 
on the fraction of halo brown dwarfs derived by Graff and Freese (1996a) is not valid in 
the present context. 
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As a matter of fact, recent N-body simulations of hierarchical cluster- 
ing produce dark matter halos with a universal Navarro-Frenk-White (1997) 
(NFW) density profile. However, this profile gets modified by the baryonic 
infall (associated with the disk formation) (Blumenthal et al. 1986; Mo, Mao 
& White 1998). In the case of Milky- Way- type galaxies, over galactocentric 
distances ranging from a few kpc up to nearly 50 kpc the NFW density pro- 
file is turned into a SIS profile, which gives rise to the observed flat rotation 
curves. Manifestly, this result contradicts the Fall-Rees assumption of a SIS 
profile before baryonic infall, and so invalidates the above discussion. For 
instance, eqs. (1) and (7) imply that the fractional abundance of baryonic 
matter increases linearly with the galactocentric distance and should even- 
tually dominate the mass density, but numerical experiments disprove this 
behaviour. In addition, self-consistency of the model (pd < p ) demands that 
the radius of the dark halos should obey the constraint r < 83 oT x l 2 kpc, but 
again there is good observational evidence that dark halos are considerably 
more extended (Zaritsky et al. 1997) 21 . 

Yet, the main point of the present dark matter scenario - namely the 
survival of H 2 at large galactocentric distances, which gives rise to the dark 
cluster formation via efficient cooling and fragmentation - is independent of 
the specific halo density profile. So, the qualitative features of the model in 
question are expected to remain true even within the contemporary picture 
of galaxy formation. Unfortunately, before this issue is settled no prediction 
about the actual distribution of the dark clusters in galactic halos can be 
made. 

2.2 Properties of the dark clusters 

Let us summarize the main properties of the dark clusters, as suggested by 
the above considerations. 

Dark clusters resemble both morphologically and dynamically globular 
clusters, apart from the obvious difference of being made of brown dwarfs 
rather than ordinary stars. This fact entails in turn further, crucial differ- 
ences, which we are now going to discuss. 

21 Besides, it has become clear that the formation of globular clusters does not proceed 
in accordance with the Fall-Rees theory (see e.g. Fall & Zhang 2001; Van den Bergh 2001; 
Beasley et al. 2002). 
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A difference between globular and dark clusters concerns their mass spec- 
trum. In fact, globular cluster masses exhibit a small scatter around a pre- 
ferred value ~ 10 5 M Q . Within the Fall-Rees model, this circumstance is 
naturally traced to the formation mechanism. Indeed, the permanence of 
the corresponding bubbles for a long time in virial equilibrium results in the 
imprinting of the associated Jeans mass, which is just ~ 10 5 M when eq. 
(10) is corrected for the mass loss. On the other hand, the bubbles that 
generate the dark clusters cool monotonically, and so no preferred mass scale 
gets singled out. Accordingly, we expect the dark cluster mass spectrum to 
be much wider. However, disruptive effects - like evaporation, encounters 
(both among themselves and with globular clusters), tidal disruption and 
spiralling motion towards the galactic centre (brought about dynamical fric- 
tion) - strongly constraint their mass range 22 . Specifically, only dark clusters 
in the mass range 3 • 10 2 — 10 6 M Q are expected to survive all these effects, 
and therefore to still populate the outer galactic halos today (De Paolis et al. 
1996, 1998a). We stress that this is possible only because the dark matter 
model in question predicts that dark clusters do not populate the inner halo. 
Typical values of the dark cluster radius lie in the range 1 — 10 pc. Moreover, 
all dark clusters more massive than ~ 10 5 M Q have presumably entered the 
phase of core collapse. 

Another difference between globular and dark clusters concerns their gas 
content. Because the process of star formation is highly inefficient, at least 
60% of the original amount of gas does not get transformed into stars (Scalo 
1985). Within globular clusters, such a leftover gas is expelled by stellar 
winds and shock waves driven by supernova explosions. The case of dark 
clusters is different. Since practically no nuclear process and no evolution 
occur in the brown dwarfs, the leftover gas - which is mainly Hi - remains 
trapped inside the dark clusters in the form of self-gravitating clouds. Indeed, 
we expect these gas clouds to provide the leading mass contribution to the 
dark clusters. 

Although these clouds are primarily made of H 2 , they are expected to be 
surrounded by a layer of neutral hydrogen HI and a ionized "skin", owing to 
the interaction with the diffuse photon background (De Paolis et al. 1996, 

22 A general analysis of these constraints has been carried out by Carr and Lacey (1987), 
Moore and Silk (1995), and Carr and Sakellariadou (1998). As they assumed that the 
dark cluster distribution extends all the way down to the galactic centre, their results do 
not apply to the present model. 
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1998a). As we shall see, the presence of HI and ions in the outer region of 
the clouds plays an important role in some subsequent considerations. 

An analysis of the cloud thermal balance shows that they are very cold, 
with a central temperature ~ 10 K. Moreover, Gerhard and Silk (1996) have 
demonstrated that for realistic values of their parameters the clouds should 
survive evaporation and collisional disruption. Typical values of the cloud 
mass and radius are ~ 10 _3 M Q and ~ 10 _5 pc, respectively. 

Mainly in view of the observational implications, an important question 
concerns the present dust-to-gas ratio of the clouds. Unfortunately, a clear- 
cut answer cannot be given. It looks natural to suppose that the clouds 
formed from almost primordial gas. But even if they were essentially dust 
free at the beginning, their interaction with the interstellar matter - when 
they periodically pass through the galactic disk - should have produced a 
sizable dust-to-gas ratio today, as pointed out by Kerins, Binney and Silk 
(2002). Accordingly, it has been suggested that the clouds should be opaque 
to visible light (Gerhard & Silk 1996; Kerins, Binney & Silk 2002). However, 
this conclusion may not be true. For instance, Draine (1998) has argued 
that dust grains could have sedimented in a small core, so that a cloud 
would be essentially transparent at optical wavelenghts. In the following, 
both possibilities will be considered. 

Cloud stability is a critical issue. Although a completely satisfactory 
treatment is still lacking, various mechanisms have been envisaged that can 
stabilize the clouds against gravitational collapse 23 . A possibility inves- 
tigated by Gerhard and Silk (1996) is that the clouds get stabilized by the 
gravitational field produced by the brown dwarfs clumped into the dark clus- 
ters. Alternatively, it has been suggested that cosmic-ray heating can balance 
either molecular cooling in dust-free clouds (Sciama 2000) or cooling by dust 
emission in dusty clouds (Lawrence 2001). Finally, Wardle and Walker (1999) 
have argued that stability can be achieved via sublimation of liquid or solid 
hydrogen, which is expected to be present in clouds with mass < 0.02 M & . 

In turn, it looks remarkable that the presence of gas clouds in the dark 
clusters can affect the properties of the brown dwarfs. Indeed - as shown 
by Hansen (1999) - the cloud-brown dwarf interaction can give rise to a 

23 As we will see later, the clouds are likely to be magnetized. However, it is not clear 
whether this fact helps to stabilize them, since the energy balance depends critically on 
the configuration of the magnetic field. 
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low-entropy accretion process 24 (Lenzuni, Chernoff & Salpeter 1992). Ac- 
cordingly, the brown dwarf mass gets substantially increased - up to 0.3 M & 
- while the central temperature stays below the hydrogen-burning threshold. 
Although it would then be more appropriate to talk about beige dwarfs - as 
indeed suggested by Hansen - we will continue to use the name brown dwarfs 
for simplicity, but it should be kept in mind that their mass can be as large 
as 0.3 M Q . 

Finally, we point out that also binary brown dwarfs should be clumped 
into the dark clusters. In fact - much in the same way as it occurs for 
ordinary stars - also in this case the fragmentation process is expected to 
produce a large fraction (up to 50% in mass) of binary objects 25 . Because 
of dynamical friction on the gas clouds, the overwhelming majority of binary 
brown dwarfs become so close (hard) today that they cannot be resolved 
in current gravitational microlensing experiments towards the Magellanic 
Clouds (more about this, later) 26 (De Paolis et al. 1998b). 



3 Observational implications 

The formation mechanism discussed in Subsect. 2.1 obviously provides a 
rationale for the existence of dark clusters with the properties stated in Sub- 
sect. 2.2. Yet, a more pragmatic attitude can also be taken (Gerhard & Silk 
1996). Indeed - regardless of the specific formation process - one can simply 
suppose that dark clusters of brown dwarfs 27 and cold H 2 clouds account 
for baryonic dark matter and populate the halos of normal spiral galaxies at 
galactocentric distances > 10 — 20 kpc. 

What are then the observational signatures of this scenario? Below, we 
discuss a few most important effects. 

7- ray emission - At the qualitative level, the situation looks quite simple. 
Very high-energy (E > 1 GeV) cosmic-ray protons travelling in the halo of 
a normal spiral galaxy produce (in particular) neutral pions upon scattering 

24 We stress that this process is likely to occur in the quiet environment inside the dark 
clusters. 

25 Binary brown dwarfs tend to concentrate in the dark cluster cores owing to the mass- 
stratification instability (Spitzer 1987). 

26 We remark that the energy released during the hardening process of binary brown 
dwarfs is efficiently radiated away by the clouds. 

27 Recall that by brown dwarfs we actually mean beige dwarfs throughout the paper. 
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on the gas clouds. And of course the pions rapidly decay into photons. So, a 
7-ray emission is expected, whose flux is proportional to the colomn density 
of the clouds. 

Unfortunately, a quantitative analysis is difficult, since nobody knows 
how the cosmic-ray protons propagate inside the dark halos 28 . Hence, some 
further assumptions are definitely needed at this point. 

A strategy to proceed can be sketched as follows, focusing the attention 
on the Milky Way (De Paolis et al. 1999, 2000). First of all, recall that in 
the Milky Way disk cosmic rays are produced by stellar winds and supernova 
explosions, and tend to escape into intergalactic space. However - owing 
to the inhomogeneities of the disk magnetic field over scales 10~ 6 — 10 2 pc - 
they undergo a diffusion process, which gives rise to a temporary confinement 
controlled by the escape time from the disk. Next, observe that - within the 
present dark matter model - a similar situation is expected to occur in the 
Milky Way halo as well. For, we have seen that the gas clouds - with a photo- 
ionized "skin" - have typical size ~ 10~ 5 pc, and they are clumped into dark 
clusters having typical size ~ 10 pc. Thus, we expect inhomogeneities of the 
halo magnetic field of the same kind as the ones occurring in the disk 29 . 
This circumstance allows for an estimate of the cosmic-ray escape time from 
the halo. It turns out that cosmic-ray protons with energy E < 10 3 GeV 
have an escape time larger than the Hubble time. Still, their spectrum goes 
like E~ a (2 < a < 3), and so just these protons give the leading contribution 
to the cosmic-ray proton flux in the halo. Therefore, most of the cosmic- 
ray protons produced in the Milky Way disk should still be trapped inside 
the Galactic halo today. As a result, we can estimate the average density 
of cosmic-ray protons in the halo, which turns out to be 30 ~ 0.1 eVcm~ 3 , 

28 We stress that - contrary to the practice used in the cosmic-ray community - by halo 
we mean the almost spherical galactic component which extends well beyond ~ 10 kpc, 
and not the thick disk. 

29 This picture is supported by Kronberg's (1994) suggested existence of a cosmic back- 
ground magnetic field with strength ~ 1 fj,G, which is motivated by the fact that magnetic 
fields with this strength are found nearly everywhere, regardless of the actual density and 
composition of the corresponding region. Moreover, the inhomogeneities of the magnetic 
field inside a dark cluster are expected to be produced both by the ionized envelop of the 
clouds and by the turbolent dynamo process likely present in the brown dwarf coronae 
(Drake et al. 1996). 

30 This value is consistent with the upper bound derived from EGRET observations 
(Sreekumar et al. 1993). 
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namely roughly one-tenth of the disk value. Such a high cosmic-ray proton 
flux gives rise to a potentially detectable r y-ray flux from the Milky Way halo 
through the above-mentioned mechanism. Owing to the very poor angular 
resolution of present-day 7-ray detectors, this flux is indistinguishable from a 
truly diffuse emission from the Galactic halo. Although its intensity depends 
on various somewhat uncertain parameters, an order-of-magnitude estimate 
yields for the integrated flux above 1 GeV (De Paolis et al. 1995a,b) 

$ 7 (> IGeV) = 1(T 7 - KT 6 7 cm" 2 s~ l sr' 1 . (12) 

We stress that an independent calculation (Kalberla, Shchekinov & Dettmar 

1999) - based on different assumptions - leads to the same conclusion. 

In spite of the fact that the flux in eq. (12) is slightly less than the 
observed diffuse extragalactic flux (Sreekumar et al. 1998), a wavelet-based 
statistical analysis can discriminate between a Milky Way halo emission and 
an extragalactic flux. In 1998 such an analysis has been carried out for 
EGRET data and has led to the discovery of a diffuse 7-ray emission from 
the Galactic halo (Dixon et al. 1998). Remarkably enough, the observed 
flux is in agreement with both the intensity and the spatial distribution of 
the emission predicted by the considered dark matter model, provided that 
a moderate halo flattening is allowed 31 (De Paolis et al. 1999, 2000). 

It goes without saying that the future planned satellite missions AGILE 
and GLAST will play a crucial role in settling this issue. 

As is well known, the Milky Way is a typical normal spiral galaxy, and so 
we expect such a 7-ray halo emission from nearly all normal spiral galaxies. 
In particular, observing the 7-ray flux from the halo of Andromeda galaxy 
will be a challange for the next generation 7-ray detectors (De Paolis et al. 

2000) . 

A final comment is in order. When high-energy cosmic-ray protons scatter 
on the gas clouds, also charged pions are produced, which ultimately give rise 

31 As already pointed out, the dark matter model in question does not predict (in its 
present form) how dark clusters are distributed in the Galactic halo. While it looks natural 
to suppose that their distribution follows that of the nonbaryonic dark matter - namely a 
SIS density profile - some flattening can be expected (Samurovic, Cirkovic & Milosevic- 
Zdjelar 1999). In the present calculation, the baryonic halo is modelled as a flattened 
spheroid — which becomes a SIS in the limit of spherical symmetry - with the flattening 
parameter left free. Because nonbaryonic dark matter dominates over the baryons, a 
moderately flattened baryonic dark halo is consistent with the recent evidence of a nearly 
spherical Galactic halo (Ibata et al. 2001). 
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to a flux of high-energy electrons 32 33 . Were these electrons to escape from 
a cloud, two further effects would come about and should be addressed. 
One is emission of synchrotron radiation in the halo magnetic field. The 
other is inverse Compton scattering against CMBR photons, leading to a 
soft X-ray flux. Although in either case the intensity is expected to be 
subdominant 34 , the latter effect might be disentangled from the background 
by its characteristic angular distribution. However, the relativistic electrons 
hardly escape from a cloud - owing to its ionized outer envelope - and so 
these effects fail to provide a signature of the dark matter scenario under 
consideration. 

X-ray emission - As already emphasized, brown dwarfs are expected to 
possess a coronal soft X-ray emission. More explicitly, Kayshap et al. (1994) 
have quantified this flux in ~ 10 27 erg s 1 in the 0.1 — 10 keV energy range. 
Consequently, the possibility arises to discover the dark clusters in X-ray 
searches. 

A thorough analysis has shown that this should in fact be the case for dark 
clusters with mass as large as 10 5 M , provided that the fractional abundance 
of brown dwarfs in the dark clusters is sizable (De Paolis et al. 1998c). 
Unfortunately, it is impossible to tell a priori whether this circumstance is 
indeed realized, and so an intrinsic uncertainty affects the present discussion. 
More specifically, dark clusters with mass ~ 10 5 M & can contribute to the 
new population of faint X-ray sources advocated by Hasinger et al. (1993) 
and by McHardy et al. (1997), whereas dark clusters with mass ~ 10 6 M can 
be observed as resolved sources with the future planned satellite missions. 

We stress that these results were obtained in 1998, when only ROSAT 
data were available. The impact of the Chandra and XMM-Newton missions 
on this issue has still to be investigated. 

Infrared emission - A different kind of dark matter search addresses the 

32 Positrons disappear because of annihilation into photons, thereby increasing the 
previously-considered 7-ray flux. 

33 Also high-energy neutrinos arise in this process, but a preliminary analysis (Ingelman 
& Thunman 1996) has shown that this flux docs not compete with that of neutrinos 
produced in the Milky Way disk. 

34 The 7-ray flux discovered by Dixon et al. (1998) is slightly smaller than the diffuse 
extragalactic flux, which is in turn considerably smaller than the flux from the Milky Way 
disk (Hunter et al. 1997). Because the branching ratio is obviously independent of the 
environment, also the electron flux from the halo should be considerably smaller than the 
corresponding electron flux from the Milky Way disk. 
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infrared emission from very low-mass stars in the dark halos of various galax- 
ies, including the Milky Way. Although strong constrains have been set on 
the abundance of red dwarfs, dark halos dominated by old brown dwarfs are 
still a viable possibility (Boughn & Uson 1995; Graff & Freese 1996b; Gilmore 
& Unavane 1998). In this respect, it should be stressed that all the analyses 
performed so far rest upon the strong assumption of a smooth distribution of 
low-mass stars within a canonical SIS halo model 35 . Yet, we have seen that 
within the present dark matter scenario things are different. Brown dwarfs 
are clumped into the dark clusters, and this very fact reduces the expected 
infrared flux (Kerins 1997a,b). In addition, we expect gas clouds - rather 
than brown dwarfs - to give the leading mass contribution to the dark clus- 
ters, and so the resulting infrared flux gets further reduced 36 . Therefore only 
future observations might detect the infrared emission in question. 

Infrared emission from both dust-less (Sciama 2000) and dusty (Lawrence 
2001) clouds has been considered in a different context, with the conclusion 
that only for clouds close to the Milky Way disk can the resulting flux lie 
above the threshold of the infrared detector SCUBA. 

Occultation effects - As discussed in Subsect. 2.2, there are good rea- 
sons to expect that the gas clouds clumped into the dark clusters are opaque 
at optical wavelenghts. Consequently, they can be detected by looking for 
occultations of background stars (Gerhard & Silk 1996). Quite recently, 
Kerins, Binney and Silk (2002) have suggested that the data sets of gravita- 
tional microlensing experiments towards the Magellanic Clouds (see below), 
the Galactic bulge and the Andromeda galaxy can also be used to search 
for occultation signatures by gas clouds. Actually, they have demonstrated 
that - for cloud parameters typical of the considered dark matter model - 
thousands of transit events should already exist within microlensing survey 
data sets. 

Optical lensing - It might nevertheless happen that the gas clouds are 
effectively transparent at optical light. In this case, the light of a background 
star can be magnified - in a symmetric fashion - when a cloud crosses its line 
of sight. The resulting light curves resemble those of gravitational microlens- 
ing (see below), apart from the fact that red light suffers a stronger magnifi- 

35 Hcnce, low-mass stars are spread out all the way down to the galactic centre with a 
r~ 2 density profile. 

36 In the case of the Milky Way yet another reduction of the infrared flux is due to the 
fact that dark clusters populate only the outer halo. 
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cation than blu light (because of Rayleigh scattering). The phenomenology 
of the corresponding chromatic events has been studied by Draine (1998) 
and by Rafikov and Draine (2001), who have computed the event rates in 
the case of background stars in the Large Magellanic Cloud (LMC). They 
also suggested that searches for gravitational microlensing could be used to 
detect optical lensing events, provided that the achromaticity constraint im- 
posed so far is relaxed (a study of the absorption features caused by the gas 
clouds would help to identify the desired events). 

Gravitational microlensing - This effect 37 towards the Magellanic Clouds 
was proposed by Paczynski (1986) as a tool to discover compact dark matter 
objects - named MACHOs (Massive Astrophysical Compact Halo Objects) 
- presumably lurking in the halo of the Milky Way. Basically, a background 
star gets magnified - in a symmetric and achromatic fashion - when a MA- 
CHO crosses its line of sight. Since 1993, the MACHO collaboration has 
detected 13 — 17 events towards the LMC, the observed optical depth being 
T i) S ~ 1.2 • 10~ 7 (Alcock et al. 2000). Only 4 events have been found by the 
EROS2 collaboration (Milsztajn & Lasserre 2001) 38 . In the following, we 
shall focus our attention on the MACHO data. 

Because a 100 % MACHO canonical SIS halo model 39 predicts T pred ~ 
5 • 10~ 7 , one would conclude that roughly 20 % of the halo dark matter should 
be in the form of MACHOs. As the theory of galaxy formation requires most 
of the galactic dark matter to be nonbaryonic, this result looks reasonable, 
but we emphasize that it relies upon two strong assumptions: all microlensing 
events are due to MACHOs, and a canonical SIS MACHO distribution. As a 
matter of fact, the former assumption is grossly violated. Besides faint stars 
in the various components of the Milky Way (thin disk, thick disk, spheroid 
and halo), also faint stars insides the LMC itself can produce microlensing 
events 40 (the latter phenomenon being referred to as self-lensing), but - as 

37 A thorough account of microlensing can be found in Mollerach and Roulet (2002). 

38 A 

comparison of the results of the MACHO and EROS 2 callaborations is not straight- 
forward, since they look at different fields on the LMC, monitoring a different number of 
stars with different observation times. 

39 Accordingly, MACHOs are smoothly spread out all the way down to the Galactic 
centre with a r~ 2 density profile. However, it should be stressed that this model is purely 
academic, since it both assumes a 100 % efficiency in the MACHO formation process and 
neglects nonbaryonic dark matter. 

40 An updated discussion of the various contributions is given by Jetzer, Mancini and 
Scarpetta (2002). 
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a rule - observations provide no information about the location of the lenses 
(unless they happen to be binary objects). Although various attempts at 
estimating the optical depth for self-lensing r se /j have been made, discordant 
results emerged (due to the poor knowledge of the mass distribution in the 
LMC), ranging from r self = (0.7 - 7.8) • KT 8 (Gyuk, Dalai & Griest 2000) 
up to r self = (0.7 - 1.8) ■ 10" 7 (Evans & Kerins 2000) 41 . Evidently no sharp 
conclusion can be drawn about the fractional abundance of MACHOs in the 
Milky Way halo, but it seems fair to state that a MACHO contribution to the 
optical depth Tmacho ~ 10~ 8 is quite consistent with present microlensing 
data. 

A set of microlensing events is also characterized by the average lens 
mass M, which is related to the event duration and depends strongly on 
the assumed galactic model. The MACHO collaboration finds M ~ 0.7 M Q 
for a canonical SIS halo model, whereas the extreme cases of a maximal 
and minimal halo yield M ~ 0.9 M and M ~ 0.5 M Q , respectively (with 
large uncertainties). Unfortunately, M does not possess any clear physical 
meaning, given that at least 5 different lens populations contribute to the 
optical depth (as already emphasized). Only by a detailed modelling of 
every population can the average mass corresponding to each population be 
estimated. According to Jetzer, Mancini and Scarpetta (2002), the MACHO 
average mass turns out to be M macho — 0.3 M & for a canonical SIS halo 
model. Obviously - within such a model - ordinary brown dwarfs are ruled 
out as a canditate for MACHOs. 

Let us now discuss how the observed microlensing events fit within the 
considered dark matter scenario. It goes without saying that - while not 
logically compelling - here MACHOs are naturally identified with brown 
dwarfs, and this fact raises two questions which will be addressed separately. 

• A question arises because then MACHOs make up a nonnegligible frac- 
tion of the halo dark matter, and so the predicted MACHO contribution 
to the optical depth t' macho might well come out too large (i. e. larger 
than - 10~ 8 ). 

41 In spite of the fact that the large value T se if — 1-4 • 1CP 7 estimated by Weinberg 
(2000) has been criticized by Gyuk, Dalai and Griest (2000), similar large values have 
been obtained on different grounds by Evans and Kerins (2000) and Zhao and Evans 
(2000). Hopefully, the new data of Van der Marel et al. (2002) on the LMC will clarify 
the situation. 
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In order to clarify this issue, we compare the dark matter model in 
question with a realistic canonical SIS halo model. Because leftover gas 
as well as nonbaryonic dark matter have to be taken into account, the 
MACHO fractional abundance in a realistic canonical SIS halo model 
cannot exceed, say, 20 %. Then the above considerations entail that 
the MACHO contribution to the optical depth predicted by a realistic 
canonical SIS halo model is t'Ij ACHO ~ 1.2 • 1CT 7 . So, the critical point 
is whether the morphological difference between the two models indeed 
entails that t'macho is lower that t'macho by roughly one order of 
magnitude. We believe that the answer is yes, as we are going to show. 

A natural expectation is that also the baryonic dark halo under con- 
sideration has a SIS density profile - but a noncanonical one - since 
only the outer halo is populated by dark clusters, hence by MACHOs. 
Thus, we can formally regard the present dark matter model as a re- 
alistic canonical SIS model with no MACHOs in the inner part. Con- 
sequently, t' macho is reduced with respect to t' macho and we actually 
get 42 t' macuo ~ (0.4 - 0.8) • 10~ 7 (Jetzer 2001). Notice that a moder- 
ate flattening of the baryonic dark halo - as suggested by the observed 
7-ray emission - would further lower t' macho . As a matter of fact, we 
expect the MACHO fraction to be less than 20 %, because the dark 
clusters should mainly consist of H 2 clouds - instead of brown dwarfs - 
and this circumstance causes an additional reduction of t' macho . Gas 
clouds are not compact enough to produce microlensing events (Hen- 
riksen & Widrow 1995), but - depending on their dust-to-gas ratio - 
can either obscurate or magnify a background star. As pointed out 
above, in the latter case the effect is chromatic. Therefore - owing to 
the achromaticity constraint imposed in current microlensing searches 
- the resulting events would simply be discarded. Now, within this 
context MACHOs are associated with the gas clouds, and so some 
would-be genuine microlensing events - produced by MACHOs - can 
either become chromatic or disappear altogheter owing to an interven- 
ing gas cloud 43 . Because the resulting events would not be observed 



42 The uncertainty depends on the minimal galactocentric distance of the dark clusters 
r = 10 - 20 kpc. 

43 A quantitative analysis of these chromatic microlensing events has been performed by 
Bozza et al. (2002). 
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in present-day experiments, yet another reduction of t' macho comes 
about. On the whole, there should be little doubt that the desired 
reduction can indeed be achieved, even if a precise estimate would be 
impossible. 

An alternative possibility has been suggested by Kerins and Evans 
(1998), who suppose that the initial mass function in the halo varies 
smoothly with the galactocentric distance. Observe that here the bary- 
onic halo model necessarily differs from the SIS. They have shown that 
in such a situation brown dwarfs fail to dominate the optical depth 
while still dominating the (baryonic) mass density. As a result, the 
risk of too large values of t' macho disappears 44 . 

• Another question concerns the predicted value of M macho, because it 
can turn out to largely exceed the brown dwarf mass, thereby preventing 
MACHOs from being brown dwarfs. 

Again, let us consider first a (noncanonical) SIS halo model. An ex- 
plicit calculation shows that the lack of brown dwarfs in the inner halo 
produces a reduction of M macho, an d we get M macho — 0.2 M (Jet- 
zer 2001). As before, a moderate flattening of the baryonic dark halo 
further lowers M macho- Yet, we have seen that within the considered 
dark matter model brown dwarfs are to be replaced by beige dwarfs, 
with masses up to 0.3 M Q . Moreover, a substantial fraction of these 
beige dwarfs should be binary systems, which are presently so close 
(hard) that cannot be resolved in current microlensing searches (De 
Paolis et al. 1998b). As a result, the effective brown dwarf mass can 
well be as large as 0.3 M . So, no problem arises. 

Within the Kerins-Evans scenario the situation is different and depends 
on the quantitative details of the model. However, M macho tends to 
come out close to the value predicted by the canonical SIS model. 
Therefore MACHOs are not brown dwarfs in this setting, and their 
nature is an oper question. A possibility could be that MACHOs are 
white dwarfs (Oppenheimer et al. 2001). Because of their low fractional 
abundance, the constraints discussed in Sect. 1 might be circumvented. 

44 Because it is not clear at this stage whether MACHOs are really brown dwarfs in this 
context (see below) , here t' macho actually denotes the brown dwarf contribution to the 
optical depth. 
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In conclusion, the considered dark matter model is fully consistent with 
present-day microlensing data 45 . 

It goes without saying that detection of chromatic microlensing events in 
future experiments would provide evidence that MACHOs are indeed sur- 
rounded by gas clouds. 

CMBR anisotropics - Dark clusters can also be discovered by performing 
high-precision measurements of the CMBR anisotropy. This strategy involves 
a sort of kinematic Sunyaev-Zeldovich effect (Sunyaev & Zeldovich 1980) 
in the gas clouds clumped into the dark clusters and can be illustrated as 
follows 46 (De Paolis et al. 1995c). Assuming for simplicity that dust effects 
can be neglected, absorption and emission processes merely involve molecular 
roto-vibrational transition lines, because of the very low cloud temperature. 
So, photons of the CMBR can be absorbed and re-emitted by the clouds. 
Since the latter photons are evidently Doppler-shifted - owing both to the 
cloud velocity dispersion inside the cluster and to the cluster peculiar velocity 
- an anisotropy in the CMBR shows up when looking towards a dark cluster. 

Of course, in order for the effect to be sizable a sufficiently large number 
of photons must be involved, and the transition line in question has to be 
optically thick. Hence, the real question is whether any such molecular line 
falls close enough to the pick of the CMBR. Although the answer obviously 
depends on the poorly-known cloud composition, even in the limiting case of 
a primordial metallicity the first rotational transition line of LiH lies very 
close to the CMBR pick and is optically thick 47 . 

Observe that for dark clusters in the Milky Way halo the resulting CMBR 
anisotropies show up on the angular scale of ~ 1 arcminute. This happens 
to be just the typical angular scale associated with CMBR anisotropies pro- 
duced - via the Sunyaev-Zeldovich effect - by clusters of galaxies (Rephaeli 
1995; Birkenshow 1999). Consequently, the same observing strategies em- 
ployed in the latter case (for a review, see Rephaeli 2001) can be used to 
detect the dark clusters in the Milky Way halo. Obviously, only the occur- 

45 Because so few microlensing events are produced by MACHOs, it is impossible to 
decide observationally whether or not MACHOs are clustered (Maoz 1994; Metcalf & Silk 
1996). 

46 We are neglecting here the kinematic Sunyaev-Zeldovich effect arising from the Comp- 
ton scattering of CMBR photons on the ionized "skin" of the clouds. 

47 This is basically due to the fact that the very low colomn density is compensated by 
the huge resonance cross-section. 
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rence of a microlensing event can tell us the actual position of a dark cluster 
on the sky - in real time, if microlensing data are analyzed on-line - and so 
a coordinated effort is required. 

A nonconventional procedure to detect dark clusters in the Andromeda 
halo - which operates on the angular scale of ~ ldegree - has been proposed 
by De Paolis et al. (1995c). 

Extreme Scattering Events (ESEs) - They are dramatic flux changes oc- 
curring - over several weeks to months - during radio flux monitoring of 
some quasars (Fiedler et al. 1987). It is generally agreed that ESEs are not 
intrinsic variations, but rather apparent flux changes caused by radio wave 
refraction when a cloud crosses the line of sight to a quasar. Evidently, in 
order to produce an ESE the cloud has to be ionized, and the radio signal fea- 
tures demand that the cloud radius and electron density should be ~ 10~ 5 pc 
and ~ 10 3 cm -3 , respectively. Assuming full ionization may look natural at 
first sight, but then the resulting electron pressure turns out to exceed that 
of the interstellar medium by a factor ~ 10 3 , thereby leading to complete 
evaporation of the cloud within ~ 1 yr. 

In 1998, Walker and Wardle (1998) claimed that the first consistent expla- 
nation of the ESEs involves an unclustered population of cold (T ~ 10 K) H 2 
clouds with radius ~ 10~ 5 pc and mass ~ 10 -3 M e , distributed in a canonical 
SIS halo and comprising most of its mass. Basically, radio way refraction is 
caused by the photo-ionized "skin", whereas the inner neutral region keeps 
the electron pressure sufficiently small. 

It looks intriguing that these clouds are practically identical to those 
clumped into the dark clusters, and so it would tempting to imagine that 
the present dark matter scenario could explain the ESEs. However, the 
cloud spatial distribution is very different in the two models and the effect 
of clustering on the ESEs has to be investigated before any conclusion can 
be drawn. 

Lyman-a clouds - A typical line of sight to a distant quasar contains a 
huge number of absorption features, brought about by intervening gas clouds 
with colomn density in the range 3 • 10 12 — 2 • 10 20 cvrT 2 . In particular, a forest 
of Lyman-a absorption lines of HI is associated with clouds having colomn 
density in the range 3 • 10 12 — 3 • 10 15 cmr 2 . Because it is difficult to es- 
tablish whether the Lyman-a lines are produced by transient filamentary 
intergalactic gas structures or by gas clouds clustered around galaxies, the 
physical nature of Lyman-a absorption systems is still controversial. While 
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both options are likely to be - at least partially - correct, Lanzetta and col- 
laborators claim that low-redshift Lyman-a lines are mostly produced by gas 
clouds inside dark galactic halos (Chen et al. 1998, 2001). Moreover, Chen, 
Prochaska and Lanzetta (2001) have estimated the cosmological equivalent 
amount of HI in these Lyman-a clouds to be &Hi,g — (2 — 8) • 10~ 3 . 

Naturally, the question arises as to whether the clouds clumped into the 
dark clusters can be identified with the low-redshift Lyman-a clouds. Su- 
perficially, the answer seems to lie in the negative. For, we have seen that 
typical cloud parameters are M ~ 1O~ 3 M and R ~ 10~ 5 pc, resulting in 
an average colomn density of ~ 10 25 cm -2 , namely ten orders of magnitude 
larger than allowed. However, such a colomn density actually pertains to H2, 
not to HI. In fact, two points should be emphasized. 

• Given that HI is present only in the outer region of the clouds 48 , 
its volume density is presumably much smaller than the average cloud 
density. 

• The reported typical values of the Lyman-a cloud colomn density are 
inferred under the implicit assumption of full HI cloud composition. 
However, this is presently not the case, and what matters here is the 
thickness of the external HI layer (rather than the cloud radius), which 
is expected to be several orders of magnitude smaller than the cloud 
size. 

Altogether, the actual HI colomn density can well lie in the allowed range. 
Furthermore, the fractional abundance of HI appears to fit within the above 
estimate if the clouds in questions indeed provide the main contribution to 
the dark baryon budget at low-redshift. Needless to say, only a detailed 
analysis of the cloud distribution predicted by the considered dark matter 
model can show whether the observed properties of low-redshift Lyman-o; 
forest lines are correctly reproduced. 

H 2 absorption lines - No doubt, in principle the most straigthforward way 
to discover the clouds clumped into the dark clusters is based upon detection 
of ultraviolet Lyman and Werner absorption lines of H2 in the quasar spectra. 
This strategy is conceptually identical to the one discussed above for the 

48 Because of this fact, the effect of dust should be irrelevant here, since dust is expected 
to dominate the inner part of a cloud. 
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Lyman-a lines of HI. However, in practice things are different. For, while 
positive detection would provide unambiguous evidence, sizable dust effects 
in the clouds would totally upset this method, and so no conclusion can be 
drawn from lack of detection. 

As a matter of fact, H 2 has already been detected in this way (Folz 
et al. f988; Ge & Bechtold 1997), but these observations suffer from a 
severe confusion problem. Nevertheless, the FUSE satellite (Sembach et al. 
2002) should collect high-quality data on the H 2 absorption lines, and so a 
clarification of this issue is to be expected. 

Clusters of galaxies - As is well known, clusters of galaxies contain a 
large amount of hot X-ray emitting gas. Besides thermal Bremsstrahlung, 
the intracluster gas also produces heavy-element recombination lines. This 
fact implies that at least part of the gas must have been processed inside 
the galaxies and subsequently ejected. Actually, specific models of galactic 
chemical evolution correctly account for the observed metallicity of the intr- 
acluster gas, but its predicted total amount invariably turns out to be about 
one order of magnitude lower than observed (Matteucci & Vettotani 1988; 
Ciotti et al. 1991; Metzler & Evrard 1994). So, it appears that roughly 90 % 
of the intracluster gas cannot be galactic in origin. 

To the extent to which the present dark matter scenario describes the 
cluster galaxies, such a conclusion is not necessarily true. In fact, what 
the above galactic evolution models implicitly assume is that only luminous 
baryons are present. However, a nontrivial amount of dark gaseous baryons 
lurking in galactic halos provides an additional supply, which can be pro- 
gressively transferred to the intracluster medium via ram-pressure stripping 
and galaxy-galaxy interaction 49 . Moreover, an order-of-magnitude estimate 
shows that a large fraction of the intracluster gas can be accounted for in 
this way. Accordingly - following a previous analysis 50 by David (1997) - 
one can conclude that the fractional baryonic content of individual galaxies, 
groups and rich clusters is nearly a constant (independent of scale). 

49 The relevance of galactic outflows for clusters of galaxies has been emphasized by 
Binney, Gerhard and Silk (2001). 

50 This analysis was based on the preliminary evidence - then turned out to be wrong - 
that MACHOs make up - 50% of the Milky Way halo. 
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4 Conclusions 



After some introductory remarks about the observational evidence for dark 
matter, we have shown that the need for baryonic dark matter is today even 
more solid than it was in the past. What is still missing, however, is some 
clear-cut observational evidence about the specific form taken by the dark 
baryons. 

Nevertheless, recent developments in computer simulations of galaxy for- 
mation as well as in the understanding of low-redshift Lyman-a clouds have 
provided valuable hints. When they are combined with metallicity and back- 
ground light constraints, the list of baryonic dark matter candidates in galax- 
ies gets dramatically shortened. 

Dark clusters of brown dwarfs and cold H 2 clouds in galactic halos cer- 
tainly look as a natural possibility. We have discussed how they are expected 
to form, the properties they presumably should have and the observations in 
which they are likely to show up. 

A good deal of theoretical and observational work has still to be done 
before claiming that the present dark matter scenario is correct. But - we 
believe - the stakes are worth the effort. 
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